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Abstract 

> , 

CN , We quantify the level of polarization of the atmosphere due to Zeeman splitting of oxygen in the 

' Earth's magnetic held and compare it to the level of polarization expected from the polarization 

■ of the cosmic microwave background radiation. The analysis focuses on the effect at mid-latitudes 

, and at large angular scales. We find that from stratospheric balloon borne platforms and for 



o 



observations near 100 GHz the atmospheric linear and circular polarized intensities is about 10" 
and 100 x 10~^ K, respectively, making the atmosphere a negligible source of foreground. From the 
'q^i ground the linear and circular polarized intensities are about 10~^ and 100 x 10~^ K, making the 

^ I atmosphere a potential source of foreground for the CMB E (B) mode signal if there is even a 1% 

i-( ' (0.01%) conversion of circular to linear polarization in the instrument. 

1 Introduction 

^ I Ground and balloon borne experiments are attempting to detect and characterize the cosmic mi- 

' crowave background (CMB) polarization signal. These efforts will continue into the foreseeable 

future and their success depends on a thorough understanding of the sources of foreground emis- 
sion. The atmosphere is a source of foreground emission because the magnetic field of Earth causes 
a Zeeman-splitting of the energy levels of oxygen and emission due to transitions between these 
levels occurs at the frequency band where the CMB is most intense. The Zeeman transitions are 
100% polarized and it is therefore necessary to quantify to what extent polarized emission and 
absorption in the atmosphere affects CMB polarization experiments. 

Tinkham and Strandberg [1, 2] and Townes and Schawlow [3] have derived the energy levels 
of the ground state of the oxygen molecule with and without the presence of a magnetic field, 
Lenoir [4], and Rosenkranz and Staelin [5] have discussed polarized emission of oxygen in the 
Earth's magnetic field. Keating et al. [6] found that at a frequency of 30 GHz the fraction of linear 
polarization due to this emission should be less than 10 x 10^ K. In this paper we give concrete 
predictions (rather than an upper limit) for the level of polarization due to the atmosphere over a 
broad range of frequencies and compare this level to the signals expected from the E and B modes 
of the CMB polarization. 
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Figure 1: The zenith brightness temperature of the US standard atmosphere [7] at an altitude 
of 5 km and relative humidity of 10% (top panel) and the spectrum of the fluctuations in the 
cosmic microwave background radiation for an RMS fluctuation of 50 /jK (bottom panel), both as 
a function of frequency. Atmospheric windows near 30, 90 and 150 GHz make observations of the 
CMB attractive at these frequency bands. 

2 Atmospheric Polarization 

Figure 1 shows the spectrum of the fluctuations of the CMB near its peak and the zenith effective 
temperature of the US Standard Atmosphere [7] at a 5 km high site. There are strong oxygen lines 
at ~ 60 and 118.8 GHz, a water line near 180 and 325 GHz, and a number of lower intensity ozone 
lines. Of these three molecules, oxygen is the only molecule that exhibits Zeeman splitting and has 
polarized emission. 

At its electronic ground state the oxygen molecule has 12 valence electrons, as shown in Fig. 2. 
The two electrons at the highest energy level pair up with their spins parallel, while the spins of all 
other electrons combine to give spin zero. Therefore the total spin quantum number of the molecule 
is iS = 1. At the mm- wave frequency band the transitions of the oxygen molecule are between 
rotational states of the electronic and vibrational ground state. The rotational angular momentum, 
with a quantum number A^, couples with the spin to give a total angular momentum state J. The 
resulting energy levels and rules for the allowed transitions between them are summarized by Lenoir 
[4]. Figure 2 shows the transitions between the first two rotational levels, = 1 and 3. The 118.8 
GHz line is the transition between a J = and a J = 1 states, both of them with A^ = 1. There 
are a number of transitions with frequencies around 60 GHz, which give the strong line between 55 
and 65 GHz shown in Fig. 1. 
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Figure 2: (a) The energy and spin configuration of the 12 valence electrons of the oxygen molecule 
in its ground state. The two highest energy electrons, which are in the 2iTg orbital, have their spins 
parallel giving a net spin S" = 1 to the molecule. The lower energy orbitals, which are in order 
the IvTg, 3a g, 2au, and lc7„, and the 4 inner shell electrons (not shown) give a net spin of zero, 
(b) The first two rotational energy levels, N=l and 3, of the oxygen molecule in the electronic 
and vibrational ground state. The N and J pairs give the rotational and total angular momentum 
quantum numbers. The numbers near the arrows give the frequency of the transition in GHz. The 
Zccnian splitting is not shown. 



In the presence of the Earth's magnetic field the 2 J + 1 degeneracy of each of the J states is 
broken giving rise to a Zeeman splitting of the energy levels and therefore to many new transitions. 
Expressions for the magnitude of the energy split as a function of the quantum numbers and the 
magnetic field strength and for the selection rules for transitions between the states are given by 
Lenoir [4]. (Notice an inconsistency in the notation of the energy splitting between Lenoir and 
Townes and Schawlow [3]; we follow the notation of Towncs and Schawlow, which is the more 
prevalent in the literature.) For the typical strength of the Earth magnetic field the frequency shift 
of the transitions due to the Zeeman splitting is on the order of 1 MHz. 

For simplicity we will focus our attention on the 118.8 GHz line, although our analysis is valid 
for all of the oxygen lines. In the presence of a magnetic field the {N,J) = (1,1) state is split 
to three states an M = — 1, 0, and +1 states ordered in increasing energy, as shown in Fig. 3. 
The transition with AM = is called the tt transition and the transitions with AM = ±1, where 
AM = Miow energy ~ ^high energy: are Called a±. These transitions are 100% polarized and their 
intensity, and the orientation of the polarization is determined by the angle between the line of 
sight and the magnetic field direction. When the line of sight is perpendicular to the magnetic field 
the two a lines are linearly polarized both in the same direction and parallel to the direction of 
the magnetic field, and the tt line is linearly polarized perpendicular to the field. When the line 
of sight is parallel to the magnetic field the tt line completely disappears and the two a lines are 
circularly polarized in opposing orientations. A more detailed discussion of intermediate cases is 
given by Lenoir [4]. Radiation propagating through the atmosphere at frequencies near 120 GHz 
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(or near any other oxygen line) will become polarized because of differential polarized absorption 
and emission by the a and tt lines. 



2.1 Order of Magnitude Estimate 

It is not difficult to give an order of magnitude estimate of the polarized intensity of the 118 GHz 
oxygen line in emission due to the Zeeman splitting. Since we are mainly interested in frequencies 
at which the atmosphere is not opaque, we can neglect Doppler broadening. We further assume 
that the line profile is collisional with no mixing, that is Lorcnzian. For the 118 GHz line this 
assumption is valid because the relative strength of the three Zeeman components changes both 
with polarization and viewing geometry. Any mixing between two transitions is reciprocal in the 
sense that it acts oppositely on them, and isolation of one transition would imply a non-physical line 
shape, such as negative absorption. One Zeeman component can not couple to another component 
that might disappear to a different observer. It can only couple to another component of the same 
AM. (See also the discussion following Eq. 3 of Rosenkranz and Staelin [5].) Hence the emission 
from the three components is incoherent and we simply sum the Stokes vectors that are associated 
with each of them 



(1) 



where S is the appropriate Stokes vector, v is the frequency of interest, 9 is the angle between the 
line of sight and the direction of the magnetic field, and /x is the relative strength of the transition, 
which is proportional to the square of the matrix element between the energy states. The Lorentzian 
functions are given by 



L{u) = Jo 
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with 



for TT 

vq — ^ for (T+ , (3) 
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where uq is the frequency of the AM = transition, S is the magnitude of the shift in frequency, 
and Iq is a normalizing constant proportional to the temperature of the atmosphere. The Stokes 
vectors of the lines are 
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and {fi7T,fJ-a+,fJ-a^) = (2, 1, 1). 

We estimate the magnitude of the effect for two limiting cases of observation parallel to and 
perpendicular to the direction of the magnetic field. In the limiting case of an observation parallel 
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Figure 3: Zeeman splitting of the {N,J) = (1,1) energy level to three levels in the presence of a 
magnetic field. The AM = transition is called tt and the AM = ±1 transitions are called a±, 
where AM = M - M'. 



to the field, for which ^ = 0, we have 



( 2 \ 




V±2; 
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and the null vector for S-,^, which gives two circularly polarized lines from the a components and 
no TT component. For the case of observation perpendicular to the field, for which 9 = 7r/2, we have 



ScF± — 



and 
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(8) 



which gives orthogonal linear polarizations for the u and tt components. 

In the limit where 7, (5 ^ {v — vq), which is always the case for altitudes of less than 50 Km 
above Earth, the polarized intensity for ^ = 7r/2 is purely a Q Stokes state and is given by 



/>,e = 7r/2) = g(z/)~i/t(-4 



(9) 



where It is the total intensity (the I Stokes parameter). The polarized intensity is quadratic in 
the ratio of 5 to the distance in frequency from line center v — vq. for u = 140 GHz we have 
V — vq = 2{) GHz and with 6 = 1 MHz, corresponding to a field of ~ 1 Gauss, and an atmosphere 
brightness temperature of 10 Kelvin we find /p ~ 10 x 10~^ Kelvin, similar to the estimate of 
Keating et al. [6]. 
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For the geometry where ^ = the polarized intensity is purely a V Stokes state, i.e. it is fully 
circular, but with a much larger amplitude 

Ip{u, e = 0) = Viu) ~ Alt (—^) ■ (10) 

The polarized intensity is now linear in the ratio of 5 to the distance in frequency from line center, 
for u = 140 GHz we have Ip ~ 10~^ Kelvin. Even a 0.1% conversion of circular polarization to 
linear in the instrument will give rise to signals that are significant compared to the cosmological 
signals. 

2.2 A More Exact Calculation 

The brightness temperature matrix of the atmosphere T^ at a frequency u is given by [4] 

Tb(z) = exp {-G(z)} Tb(0) exp {-G*(z)} + Tb(0) [I - exp {-2A(z)}] , (11) 

where G(z) is the propagation matrix at a height z above Earth, A is the attenuation matrix and 
I is the identity matrix. We use the computer code of Rosenkranz and Staelin [5] to calculate the 
polarized radiative transfer of the cosmic microwave background radiation through two standard 
atmospheres [7, 8]: The 1976 US Standard Atmosphere, which is characteristic for latitudes of 
about 40N with total water vapor of 14 precip. mm, and a winter atmosphere at a latitude of 60N 
with total water vapor of 4 precip. mm. Magnetic field information comes from the IGRF 1995 
magnetic field model [9] . Tabic 1 gives the effective temperature of each of the stokes parameters 
(in thermodynamic temperature) at three different frequencies and two altitudes of observations. 
The observation sites are at a longitude of 140 W and at either 40N or 60N, as appropriate, and 
observations are pointed toward North at an elevation of 45 degrees. The calculations include 
contributions due to the 60 GHz oxygen band, for which we calculate the line mixing using the 
parameters that were measured by Liebe et al. [10]. These measured values of line mixing (and 
line broadening) correspond to the combination of Zeeman components within each line, and hence 
are appropriate to the frequencies considered here, which lie far outside the region of Zeeman 
structure. However, all of the frequencies considered here lie within 90 GHz of the line centers, so 
duration-of-collision effects are considered to be negligible. 

3 Discussion 

The results of the calculations Tising the radiative transfer code and the order of magnitude estimate 
generally agree and show that the magnitude of the circular polarization is significantly larger than 
that of the linear polarization component. Although the CMB itself is not expected to be circularly 
polarized, a 1% conversion of circular to linear polarization in the instrument at a frequency of 90 
GHz would give rise to a 3 — 6 x 10~^ K signal, depending on the location of observation. Such 
a signal is of the same order of magnitude as the strongest 'E mode' signal and much larger than 
the strongest 'B mode' signal expected from the CMB, as shown by the bars in Figure 4. However, 
the Figure also illustrates that the spatial scale where the CMB peaks and the spatial scale that is 
relevant for the calculation in this paper arc quite distinct. (Spatial scales are quantified in terms 
of £ ~ 180/0, where 9 is the angle between two lines of sight in degrees.) Under the assumption 
that the atmosphere is well mixed, an assumption that underlies all atmospheric modeling, spatial 
variations of the atmospheric polarization arise because of variations in the magnetic field along 
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Figure 4: The angular power spectrum of the CMB E (solid) and B (dashed) mode polarization 
signals in a ACDM cosmology with tensor to scalar ratio r = 0.01. The top (bottom) short bars 
near ^ = show the level of atmospheric polarization at large angular scales that is expected for a 
1% (0.01%) conversion of circular to linear polarization in an instrument making observations from 
the ground at 90 GHz and at a latitude of 60N. 

different lines of sight. However, the angular resolution of the IGRF model provides description of 
spatial variations at ^ ~ 1 for all latitudes except close to the poles. We are currently investigating 
the magnetic field information close to the poles and at scales smaller than those provided by IGRF 
and plan to report on it in a future publication. 

The polarization of the atmosphere should be considered carefully by ground based experiments 
that attempt to measure the CMB B mode polarization signal. As Fig. 4 shows, a 0.01% conversion 
of circular to linear polarization in the instrument gives rise to an atmospheric polarization signal 
that is larger by x2 than the largest B signal expected with a cosmology that has a tensor to 
scalar ratio r = 0.01. At large angular scales the atmospheric polarization is larger than estimates 
for the polarization signal due to synchrotron [11] and Galactic dust [12]. Oblique angle reflection 
of cicularly polarized light from telescope mirrors can readily generate a 0.01% conversion to linear 
polarization. 

Because the column density of oxygen is much smaller at stratospheric balloon altitudes the 
signal observed from such a platform is smaller by about a factor of 1000 making it an unimportant 
foreground for CMB polarization experiments. 
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Figure 5: The V Stokes parameter as a function of frequency for observation through a US standard 
atmosphere from an altitude of 1 km. 



The atmospheric polarization signals that we have calculated are valid at a given frequency and 
have not been integrated over a band of frequencies; Fig. 5 shows the expected variation of the 
signal as a function of frequency between 125 and 160 GHz for an observation at an altitude of 
1 km. We do not give detailed information near the line center at 119 GHz because the line shape 
is not accurate enough at these frequencies. If the magnitude of the atmospheric polarization at 
high £ is the same as that at low £ then a strong rejection of out of band power will be necessary 
for ground based polarization experiments. 

Unlike fluctuations in the brightness temperature of the atmosphere, which are primarily due 
to water clouds at low altitudes the polarized intensity is due to oxygen and is not expected to vary 
with time, except that it may be attenuated by varying water vapor. For a given observation site 
the signal will be fixed in a particular azimuth and elevation directions, but not in right ascension 
and declination. 
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Tabic 1: Absolute value of the /, Q, and V Stokes parameters for observation through the atmo- 
sphere in Kelvin (thermodynamic units) for observation sites at two latitudes and two altitudes 
and at the given frequencies. Both sites are at a longitude of 140 W and the observation is toward 
the North at an elevation of 45 degrees. 
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